ITP-SB-02-9* 
INLO-PUB-02/9^ 

TWO-LOOP OPERATOR MATRIX ELEMENTS 
CALCULATED UP TO FINITE TERMS FOR POLARIZED 
DEEP INELASTIC LEPTON-HADRON SCATTERING 



Y. Matiounine 1 , J. Smith 1 and W.L. van Neerven 2 
1 Institute for Theoretical Physics, State University of New York at Stony Brook, Stony Brook, 

NY 11794-3840, USA 

2 Instituut-Lorentz, University of Leiden, PO Box 9506, 2300 RA Leiden, The Netherlands 

(March 1998) 



Abstract 

We present the two-loop corrected operator matrix elements contributing 
to the scale evolution of the longitudinal spin structure function g±(x, Q 2 ) 
calculated up to finite terms which survive in the limit e = iV — 4 — > 0. 
These terms are needed to renormalize the local operators up to third order 
in the strong coupling constant a s . Further the expressions for the two- loop 
corrected operator matrix elements can be inserted into one loop graphs to 
obtain a part of the third order contributions to these matrix elements. This 
work is a first step in obtaining the third order anomalous dimensions so 
that a complete next-to-next-to-leading order (NNLO) analysis of the above 
mentioned structure function can be carried out. In our calculation particular 
attention is paid to the renormalization constant which is needed to restore 
the Ward-identities violated by the HVBM prescription for the 75-matrix in 
iV-dimensional regularization. 

PACS numbers : 11.15.Bt, 12.38.-t, 12.38.Bx 
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I. INTRODUCTION 



During the last ten years there has been a lot of activity in the area of spin physics in 
particular the study of the structure functions gi(x,Q 2 ) (i = 1,2) measured in polarized 
lepton-hadron scattering (for reviews see |l| and [0]). One of the achievements was the 
determination of the anomalous dimensions of twist two local composite operators up to 
two-loop order (see [§], @|) which determine the scale (Q 2 ) evolution of the longitudinal 
spin structure function gi(x, Q 2 ). Since also the order a s contributions to the coefficient 
functions are known (see ||, ||) it is now possible to make a full next-to- leading order 
(NLO) analysis of gi(x, Q 2 ) analogous to what has been done for the spin averaged structure 
function F2(x,Q 2 ). When the statistics of the ongoing experiments improve it will be also 
necessary to investigate how the NLO predictions are modified by including yet higher 
order corrections. In particular this will be interesting for the study of the small x-region 
where gluon contributions to the anomalous dimensions and the coefficient functions are 
important. A start of this programme, which will lead to a complete next-to-next-to-leading 
order (NNLO) description of the longitudinal spin structure function, has been made in |7j 
where quark as well as gluon coefficient functions have been computed up to order a 2 . The 
three-loop anomalous dimensions are still missing. 

In this paper we want to make the first step in a three loop programme by computing the 
operator matrix elements (OME's), obtained by sandwiching the local composite operators 
between quark and gluon states, up to two-loop order including non-pole terms which are 
finite in the limit iV — > 4. Here N refers to the method of iV-dimensional regularization which 
is used to regularize the ultraviolet divergences occurring in the OME's. Using this method 
the latter divergences manifest themselves as pole terms of the type 1/(N — 4) fc . These finite 
terms are needed to renormalize the three-loop graphs contributing to the OME's. Further 
one can use the unrenormalized two-loop expressions to determine a part of the three-loop 
contributions to the OME's by inserting them into one loop diagrams. Another interesting 
feature of the calculation of the spin OME's is the appearance of the 75-matrix and the Levi- 
Civita tensor in the operator vertices. Since we use iV-dimensional regularization one has to 
find a suitable prescription to extend these objects to N dimensions. In this work we have 
chosen the HVBM prescription given first by 't Hooft and Veltman |J and worked out in 
more detail by Breitenlohner and Maison ||. One of the characteristics of this prescription 
is that the 75-matrix commutes with the other 7^ matrices for N > A. Although this is 
a consistent scheme, since it preserves the cyclicity of the traces, it violates some Ward- 
identities which would be preserved with an anti-commuting 75-matrix. To restore these 
Ward-identities one has to introduce additional renormalization factors called Z^. These 
factors will be calculated up to order a 2 for the non-singlet (r = NS) and singlet (r = S) 
operators for arbitrary spin. 

The paper is organized as follows. In section 2 we present the algebraic formulae for the 
unrenormalized OME's expressed in renormalization group coefficients. We determine the 
factors for r = NS (nonsinglet) and S (singlet) and compare our results with earlier 
calculations in the literature. The long expressions for the spin OME's, which are the results 
of our calculations, are given in Appendix A. 
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II. THE CALCULATION OF THE TWO-LOOP OPERATOR MATRIX 

ELEMENTS 



In this section we will give an outline of the calculation of the OME's up to two-loop 
order. The operators, which appear in polarized lepton-hadron scattering, can be split into 
singlet and non-singlet parts with respect to the flavour group. In leading twist (namely 
two) the non-singlet quark operator of spin n is given by 



^-\5,|ti,|ta— /*n 
U q,k 



■n-1 



S 



if){x)^D 



M2 



D^— ip(x) + trace terms 



In the singlet case there are two operators. The quark operator is represented by 

1 



1 



— I 
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n-l 



s 



ip{x)^ b ^D^ 2 ■ ■ ■ D^ n ijj(x) + trace terms 



and the gluon operator is given by 
1 



a 



•n-2 



S 



D 



,x) + trace terms 



(2.1) 



(2.2) 



(2.3) 



In the composite operators above if) and stand for the quark field and the gluon field 
tensor respectively. The in Eq. (23) represent the generators of the flavour group and 
the index a in Eq. (|2.3| ) stands for the colour. Further the above operators are irreducible 
tensors with respect to the Lorentz group so that they have to be symmetric and traceless in 
all their Lorentz indices //j. From the operators above one can derive the Feynman rules for 
the operator vertices in the standard way (see e.g. || and ||11|| ). This derivation is facilitated 
if the operators are multiplied by the source 



J, 



(2.4) 



with A 2 = in order to eliminate the trace terms in Eqs. (|2.1|) -( |2~3| ). Hence all operator 
vertices in momentum space are multiplied by a factor (A -p) n . For the computation of the 
OME's denoted by 



A 5 



(j(p)\Ot\j(p)) 



(2.5) 



with i,j = q, g we choose the Feynman gauge except for the one-loop graphs for which we 
take the general covariant gauge. For this choice the gluon propagator equals 



A£(*) 



iS, 



ab 



k 2 



(2.6) 



The matrix element (|2.5|) has to be considered as a connected Green function with the 
external legs amputated but with the external self energies of the partons j included. In this 
paper all quarks and gluons are taken to be massless and the external momentum p is off-shell 
(p 2 < 0) in order to get finite expressions for the OME's. This choice implies that the OME's 
are not gauge invariant so that they cease to be ordinary S-matrix elements. Therefore they 
acquire unphysical parts which usually vanish by virtue of the equations of motion (EOM) 
if the external legs are on-shell. Contrary to the spin averaged operators treated e.g. in 
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T2|| there is no mixing between gauge invariant (GI) or physical (PHYS) operators and 
non gauge invariant (NGI) operators (see [|l(J- Jl8|). Therefore the renormalization of the 
operators in Eqs. (|2.1|)-(fO|) is much easier than is the case for the spin averaged operators 



in 



I4 . On the other hand, due to the presence of the 75-matrix and the Levi-Civita tensors 
in the composite operators above, one has to correct for spurious terms in the OME's if one 
adopts iV-dimensional regularization. 

The calculation of the Feynman graphs corresponding to the physical operators in Eqs. 
(p.l|)-( p73|) , which are depicted in the figures in |J, proceeds in the standard way. The 
corresponding Feynman integrals reveal ultraviolet divergences which are regularized using 
the method of iV-dimensional regularization. In this way the above divergences show up in 
the form of pole terms of the type (1/s) with e = N — 4. In ||, @ it was sufficient to 
evaluate the one-loop graphs up to finite and the two-loop graphs up to single pole terms 
in order to get the second order anomalous dimensions. Here we have to include terms 
proportional to e in the one-loop expressions and the two-loop graphs have to be computed 
up to finite terms. The way to compute the two-loop Feynman integrals up to finite terms 



is presented in |]I2"|| , flIE| , and Appendix B of |19[. We used the program FORM ]2D| to do 
the necessary algebra. 

As has been mentioned above, the presence of the 75-matrix, which is essentially a four 
dimensional object, will cause some technical problems when one chooses the method of 
iV-dimensional regularization. One has to find a prescription to give a suitable definition 
valid for all space-time dimensions. We will adopt the HVBM prescription given by 't Hooft 
and Veltman || which has been worked out in more detail by Breitenlohner and Maison || . 
In order to facilitate the calculation of the OME's it is more convenient (see [^TJ, ^2|) to 
replace the term ^75 appearing in the operator vertices (see Appendix A in 0) by 

#75 = ^P- A/ W7 r - (2-7) 

Notice that this replacement is only equivalent to the HVBM prescription if a single 75- 
matrix is present in the numerator of the OME's, which is the case here. Thus only one 
Levi-Civita tensor appears in all the numerators of the Feynman integrals so that the 75 
matrix and the Levi-Civita tensor are on an equal footing. Since one has to be careful 
with the treatment of the Levi-Civita tensor, which is also a four dimensional object, we 
have to follow the procedure in Appendix B of . First one computes the numerator in the 
Feynman integral corresponding to a specific graph. This numerator contains the integration 
momenta, which are in iV-spacetime dimensions. Then one applies tensorial reduction to 
express the whole integral into tensors containing A M in Eq. ( |2.4j ) and p^ in Eq. ( |2.5| ), 
which are iV-dimensional vectors. Finally one has to project the whole matrix element on 
the tensor structure characteristic of the specific OME under consideration, which follows 
from Lorentz covariance in four dimensions. Therefore the Lorentz indices of p^ and 
with N > 4 are simply dropped. 

The advantage of the HVBM prescription is that the cyclicity of the traces is preserved. 
On the other hand it destroys the anti-commutativity of the 75-matrix. This will lead to 



a renormalization of the non-singlet axial-vector current 0^'% in Eq. (|2.1|) in spite of the 



fact that it is conserved |[23|| . Furthermore as was shown in [22] the Adler-Bardeen theorem 



| 24| is no longer true. The anti-commutativity and the Ward-identities can be restored by 
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introducing an additional renormalization constant (r = NS, S) which, when expressed 
in unrenormalized quantities, has the following form 



Zqq — 1 + CL S S £ 



(1) , 8,(1) 
11 11 



+a 2 s S 2 



-MS + <f ) - 2A>4 ,(1) 



11 



3 U ii 



(2.8) 



The hat indicates that all quantities are unrenormalized with respect to coupling constant 
a s , gauge constant £ and, in the case of the OME's, also the operator renormalization. Here 
S e denotes the spherical factor characteristic of iV-dimensional regularization 



exp 



U~1e - ln4vr 



(2.9) 



where 7^ denotes the Euler constant. Further we introduce a shorthand notation for the 
strong coupling constant so that 



Art 



9_ 

Air 



If we apply coupling constant renormalization 



l + a s S E (2(3 



(2.10) 



(2.11) 



and gauge constant renormalization 



1 + a, 



< s *( z e-) 



(2.12) 



where both are presented in the MS scheme, we obtain 



75, r 
J 11 



1 + a.S e 
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ez 
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+a 2 s S 2 £ -p z£> + z, 
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(2.13) 



In QCD (SU(N)) the renormalization group coefficients are given by 



r ll n 8 t 



^(-y-(i-0) + Wz>. 



(2-14) 



with Cf = {N 2 — 1)/2N, Ca — N, Tf = 1/2 and rif stands for the number of light flavours. 
Notice that the coefficients and z^ 2 ^ in Eqs. ( |2.8|) and ( |2 . 1 3| ) are gauge (£) independent 
in contrast to the terms proportional to e like z e qq ^ which depend on the gauge. As we 
will show later on z$ is universal but z^ 1 ' and z r '^ depend on the given quantity. In the 



non-singlet case Z qg NS can be inferred from the ratio 



jiNS.PHYS 

75,NS _ ^gg / 9 i r\ 

^gg i 5.ns.phys ' \-- > ■>) 



where ^^s.phys anc j ^5^ns,phys ^gno^g ^g physical parts (for definitions see below) of the 
OME's without and with the 75-matrix respectively. The former OME has been fully calcu- 
lated in |T2j whereas the latter OME will be presented below and in Appendix A. Notice that 



the numerator and the denominator would be identical if we applied the so-called naive (or 
anticommuting) 75-prescription. The latter method allows us to anticommute this matrix 
with the other 7^-matrices so that it reaches another vertex in the diagram and leads to 
-^gg NS = 1- Therefore the renormalization constant Z qq restores the anti-commutativity of 
the 75-matrix. The singlet coefficient which can be split into 

4 (2) = 4 S ' (2) + ^ . (2-16) 

can also be calculated by anticommuting the 75-matrix. The computation of the purely 
singlet part z qq '^ will be presented at the end of this section. 

To check our results for the sums of the Feynman graphs it is useful to have explicit 
expressions for the pole terms. Therefore we will now present the OME's expressed in 
renormalization group coefficients defined in ||. The explicit formulae for the OME's can 
be found in the Appendix A and contain terms which are finite in the limit e — > 0. Further 
it is implicitly understood that all quantities in this section, in particular the anomalous 
dimensions 7^ (i, j = q,g), are Mellin transforms even though we suppress the moment index 
n to simplify the notation. (Another way to interpret the formulae is that the OME's are 
given in parton momentum fraction (z) space when the anomalous dimensions are replaced 
by minus the corresponding Altarelli-Parisi splitting functions and the multiplications are 
replaced by convolutions.) We have written the OME's in such a way that all renormalization 
group coefficients appearing in the expressions below are given in the MS-scheme. 

In the rest of the section the superscript 5 will be suppressed because we will only discuss 
the spin OME's. Up to order a 2 a the non-singlet and the singlet OME's can be decomposed 
into 



A gg 



75^^ PHYS + 75^i;f° M ] (A • p)- 1 , (2.17) 



where A r q ^ nYS and A r q ^ OM stand for the physical and unphysical parts respectively with 
r = NS, S. The presence of the latter is due to the fact that the equations of motion are 
not satisfied. The non-singlet physical OME can now be expressed in renormalization group 
coefficients as follows 



i^< pHYs = 1 + r /2 fk N g s ' {0) + < s ' (1) - ^ 

H 1 



iqq 1 -qq "qq 



^{^(7^ (0) ) 2 -^ (0) } 
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^ U 99 ^99 ^99 V 99 99 J 

-2/3 (< s -« - + 7 S S ' (0) « S ' £ ' (1) - 



«' £ ' (1) - 49 (1) ) 



(2.18) 



The 7^ "* denote the coefficients of the order a(? +1 terms appearing in the series expansions 
for the anomalous dimensions. Using the same notation we can also express the unphysical 
part of the non-singlet OME in Eq. (|2.17|) in the aforementioned renormalization group 
coefficients 



.4 



NS.EOM 
99 



/i 



^ 2 (^) £ [7{7^ (0) C S ' (1) -2/3oC (1) 



99 



(1)7NS,(1) _ P 
Z qq U qq 



99 



The singlet OME can be decomposed as follows 

4S _ ^ NS , 2PS 
99 99 99 

where the purely singlet physical part equals 

^PS.PHYS _ fj2o2(Z^_Y\l_ (0) -,(0)1 



(2.19) 



^2.20) 



99 



+ 1/I~PS,(1) , ^(0) (1)1 , PS,(2) _ PS,(2) 
p 1 9 ^99 ^ 'qg a gq J ^ "99 Z 99 
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^ '99 99 



and the purely singlet unphysical part equals 
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PS,EOM _ -2 q2( P 
qq - " - ■' ■' 



^[i/~(o)6(i)\ + 6PS,(2) , (0)^,(1) 

^ L f \'99 S? J T 99 ^ ' 99 99 



(2-21) 



^2.22) 



7 



The next OME is 



with 



A 



1 

A • p' 



(2.23) 
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Next we need 



where 



. z (l)„(l) _ ^.(1)^(0) _ Q( ?3 



2 e 
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Finally we need 



with 
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Explicit formulae for the above unrenormalized OME's corrected up to second order, 
which contain all the finite terms independent of e, are given in Appendix A. The pole 
terms agree with the results in Eqs. (2.18)-(2.29). Since the external quark and gluon legs 
are off-shell one can directly insert the formulae in Appendix A into Feynman integrals for 
one-loop graphs. In this way one gets results for subsets of the three-loop graphs. It is 
however clear that the most difficult Feynman integrals containing non-planar diagrams, 
where all quark or gluon lines cross, remain to be done. 

The renormalized OME's are given by 



4NS _ 7 5,NS/7-l\NS JNS 
A qq - ^qq l Z ) qq A qq > 



(2.30) 



and 



Aqq — Zqq (Z ^qqAqq + (Z ^qgAgg , 



Aqg — Zgq(Z ^qgAqg + (Z ^qgAgg , 



Agq — Zg q S (Z ^gqAqq + (Z l ) ggAgq , 



Agg — Zg q S (Z X ) gq A qg + (Z r )ggA : 



gg A '-gg • 



(2.31) 
(2.32) 
(2.33) 

(2.34) 
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The operator renormalization constants are given by 



-INS 
Jqq 



(2.35) 



where 



[Z- 1 ) 



L(o) 
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NS,(1) 



2e 



Iqq 



-1\FS 
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(2.37) 



The other singlet operator renormalization constants are given by 



\Z 1 )qg — a sS £ 
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+ fl 2 9 2 [J-/I/V M°) i _(o) (o)\ _ ^(0)1 
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(2.38) 
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2s 



^gg 



(2.40) 



The above renormalization constants are chosen in such a way that the anomalous dimensions 
7j-.- are presented in the MS-scheme. Naively one would expect that the single pole terms 



in Eqs. ( [2.38| ) and ( |2.39|) are then given by — lig an d — 7oo respectively However due 



to the HVBM prescription the renormalization constants Z^- have to be modified in order 
to bring the renormalized OME's in the standard form. By this we mean that the latter 
satisfy the renormalization group equations with the anomalous dimensions 7^. This we 
will leave as a check for the reader. Notice that only the final OME's are independent of 
the prescription for the 75-matrix and the Levi-Civita tensor. 

To conclude we discuss properties of the renormalization factors Z^ r , r = NS, S in ( |2.8| ) 
computed in this paper and compare them with other results in the literature. Starting with 
the non-singlet term, which is computed via ( |2.15| ), the unrenormalized expression can be 
found in Eq. (|A11|] 



Due to the factor (— l) n (for its origin see below Eq. (|A3|) ) it can be 



split into two parts as follows 



10 



7 5,NS 

J n 



dz z 



n—l 



11 



(2.41) 



After coupling constant and gauge constant renormalization (see Eqs. ( |2.11|) , ( |2.12|) ) the 
two parts become 



75,NS,(+), 
J 11 y 



5(1 -z) + a s S £ C F [ - 8(1 - z) + s{ - 4(1 - z) ln(l - z 

-4(1-^) lnz + 2z- (1 -£)} 

+a 2 S 2 [C 2 F { - 16(1 -z)- (16 + 8z) lnz 

+ 16(1 - z) mzln(l - 2)} 

+CaC f { - -&1 - z)) - ^(1 - 2) + 8(1 - z)t{2) 



+ (-y + j^) ln z ~ 4(1 - 2) In 2 z) 
+nf C F T f {±(f(l-z)) + ^(l-z) 



+—(l-z) \nz} 



(2.42) 



and 



^ NS ' (_) (^) = a 2 s S 2 [(C 2 - ic A C F ){8(l + z)(4Li 2 (-*) + 41nzln(l + z) 

+2C(2) -ln 2 ^-31nz) -56(1-2)} 
respectively. Choosing n — 1 in Eq. Q2.41 ) we obtain 



(2.43) 



4 + e(5-(l-0 



c| {22} + c^ { -=i-f} 



^2.44) 



This is the renormalization constant which has to be multiplied with the non-singlet axial- 
vector current in Eq. ( |2.1| ) so that the latter will not be renormalized (i.e., remain a 
conserved current) if one employs the HVBM prescription. This renormalization constant 
has been calculated up to order in |22j| , [23]. Except for the term proportional to e, which 
was not presented previously, our result in Eq. Q2.44D agrees with that in the aforementioned 



literature. Recently the factor has been also computed in |^5J for the coefficient 



11 



functions corresponding to the longitudinal spin fragmentation function in e + e~-scattering. 
If we indicate the latter factor by the superscript T (timelike) and the one in Eq. ( p. 41 ) by 
S (spacelike) we have the following relation 



(2.45) 



This result demonstrates the breakdown of the Gribov-Lipatov relation p6| which is therefore 
only valid up to order a s . However for the minus component 



^s,(-), T( ,) 



^ 5 r ( - ),s ( £ ) 



(2.46) 



so the Gribov-Lipatov relation still holds at least up to second order. The above relations 
have also been reported for the spin averaged non-singlet spacelike (S) and timelike (T) 
splitting functions P^ s in Eqs. (6.37) and (6.38) of |27j where z^J is replaced by Pffl- 



Notice that there is no difference between spin and spin averaged splitting functions as far 
as the non-singlet part is concerned (see [p8~|). However there is a difference between Z^ s 
if one looks at relation (6.46) in [27]. In our case it reads 



and Pjs 



Z 5,NS,( + ),T W = Z ^S M ,S {z) + lnz f 1 dy (^S,(l) ( £ } p(0) (y) 

Z J z x V 



(2.47) 



whereas in Eq. (6.46) the In z is shifted to the integrand where it becomes In?/. Hence the 
first moments of Z^ NS ' T and Z^ NS,S are not equal anymore. It implies that the coefficient 
22 of the C 2 F part of Z^ NS ' S in Eq. ( ggj ) is changed into 12 + 16((2) for Z^ N S - T (see Eq. 
(3.46) in H). 

The calculation of Z^^ 2 \ which can only come from A^ S ' PHYS in Eq. ( 2.21[ ), proceeds 
as follows. Up to second order only one diagram contributes to this OME, which can be 
found in Fig. 3.1 of [[|. Following the same arguments as given below Eq. ( |2.15| ) one has to 
anticommute the 75-matrix to another vertex in the diagram so that Z^ s = 1. However in 
this procedure the Adler-Bell-Jackiw (ABJ) anomaly, characteristic of this graph, will also 
then be removed. In this case the expression for Eq. ( |2.21| ) becomes 



.4 



PS,PHYS,ABJ 
99 



—p 
1* 



!/I 7 (o) 7 (o)\ 

r 2 \ O '99 ! 99 J 



-{- 



Iqq 1 u qg Igq J ~ h 



99 



.(or 
99 



(0) a ( 
99 99 



'9 1 ' 



,PS,(2) 
\q 



19 99 



'99 99 



(2.48) 



Comparing the above expression with the one in Eq. ( |2.21| ) we observe that in the former 
the effective anomalous dimension has become Jgg'^ + Vqglfq with v^J = 87/(1 — z). If 
we take the first moment of this expression this effective anomalous dimension becomes zero 
meaning that the ABJ-anomaly vanishes. Therefore there is a striking analogy between 
above and zQ) in Eq. ( |2.8| ). Both constants remove anomalies i.e. the former in the case of 
the singlet current 0^ and the latter in the case of the non-singlet current 0'%. However 
the ABJ anomaly has to be restored. Hence from Eqs. Q2.21|) , (|2.26|) and Q2.48 ), we get 



12 



.4 
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PS,PHYS 

gq 



ci s S £ 



,(D 
gg 



A 



PHYS 

gg 



25 V w 



a 2 S 2 



r PS,(2) 

"gg 



(2.49) 



where the last term before the equal sign restores the ABJ anomaly in Ag g s,PHYS ' ABJ . Finally 
we obtain 



7 PS,(2) 
"<?<7 



n f C F T f 16(1 - z) + 8(3 - z) In z + 4(2 + z) In 2 z 



(2.50) 



If one takes the first moment of this expression the result becomes 3n/CV7/ which agrees 
with the one quoted in [^2 . 



Finally one could ask why it is preferable to choose the HVBM instead of the naive 75- 
prescription since in the latter case = 1. The first reason is that HVBM was originally 
proposed to obtain the ABJ-anomaly in the fermion triangle graph, which is not always 
correctly reproduced using the naive prescription (see ( |2.48|) ). The second reason is that 
the Levi-Civita tensor appearing in the OME A gq induces the HVBM prescription in the 
subgraphs containing quark loops. Therefore the naive 75-prescription does not prevent an 
additional renormalization so it is inconsistent. It is better to use a consistent procedure 
like the HVBM where all constants are fixed once and for all. 
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APPENDIX A: 



In this Appendix we present complete expressions for the two-loop corrected OME's 
computed from the Feynman diagrams depicted in ||. The second order contributions are 
calculated up to finite terms which survive in the limit e — > 0. The OME's presented here 
are unrenormalized and external self-energy corrections are included. In these expressions 
definitions of the Riemann zeta-functions ((n) and the polylogarithms Li n (z), S n ^ m (z) can 
be found in |p9| . Also the distributions (1/(1 — 2;))+ and (ln(l — z)/{ \ — z)) + and written as 
1/(1 — z) and ln(l — z)/(l — z) respectively to shorten the formulae. Note that the OME's 
given in the text are the moments of the functions listed here so 

A% = f dzz^A^z,^,-), (Al) 

where for simplicity we have not written the moment index n on the functions. Also to 
simplify the expressions we define the phase-space factor 

F = fiS £ (^V 2 . (A2) 



13 



We first split into physical and unphysical parts following the notation in Eq. ( ]2.17|) . 
The physical part is (see also Eq. ( |2.18| )) 
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/r e 
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+F C F 
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+2( 
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1 8 
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2 



1 - z 
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' 1 + ^^ + ^ + 5(l-.){-2 + C (2)}]} 



1 - z 
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1-z 
48 
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2<f(l-*)[9-16C(2)] 



+32 ln(l -z){l + z- T ^—) - 8 In z(3 + 3z - - )] 



44 44 
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1-Z' 
1 



1-2' 



- 225(1 - z)] 



31-2 

+n/CVT/[ _^(l + z _ + 85( i _ z)] } 

-z)[^-36C(2)-8C(3)] 



+ifC 2 [-56 + 12z + -^- + 5(l 
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-4 ln(l - z)(-15 + 5z + ) + 241n 2 (l - z)(l + z 



+41nz(l + 3z 
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16 lnzln(l — z) 
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1 — z' 
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--{21n(l - 2) + lnz}(l + z - — — 
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188 32 



56 
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32 4 
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+54C(3) - — C(2)1 + ln(l - z)(-78 + 622 + 
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1 — z 3 
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zY l — z 6 I — z 
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4Q 2 11Q 
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k 92 2 
xC(2)(l + 32 + 42 2 - ) + {- ln 2 (l - z) - In 3 z}(l + z - - ) 

L Z O 1 Z 

k i 92 
-2 In <(2) (3 + z + Az 2 - - ) + - In 2 2(23 + 71z + I6z 2 - - ) 

X Z X z 
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+ - In ^(38 - 642 + ) + - In z ln(l - z) (14 + 5z - - ) 
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+41nzLi 2 (l -z)(2 + Az 2 - -i— ) - 4Li 2 (l - z)(3z - — *— ) 

1 — 2 1 — 2 

-4(3 + + 42 + ^2 2 )[hi2ln(l + z) + Li 2 (-2)] 

+ (1 - 2 + 42 2 - -^)[ln 2 2ln(l -2) +21n(l -2)Li 2 (l -2) 

1 — 2 

-6Li 3 (l - 2)] + 8[ln2Li 2 (-2) - 2Li 3 (-2)](l - — *— ) 

1 — 2 

+S 12 (1 - 2)(14 - 62 + 242 2 - -^)] 
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+n / cyr / [-[2o-i9z-. ;1 
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-2(-l) n (C 2 - ^C A C F )(-{ - 8 + 8z - 41nz(l + z) 
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+8(z + — j— ) ln(l + z) [C(2) + In z ln(l + z) + 2Li 2 (-z)] 

-4(1 + z)[lnzln(l - z) + Li 2 (l - z) + Li 3 (-*)] 

18 4 

— In z(l - llz) - In 2 z(2 - 4z + -z 2 ) - 2 In 2 z ln(l + z)(3 - z ) 

3 3 1 — \~ z 

4 1 

-(12 - - + 9z - 4z 2 )[lnzln(l + z) + Li 2 (-z)] 

3 z 



-2C(3)(3 + z- 



1 + z' 



+8Si 2 (-z)(-l +3z + 



l + z' 



4(1 -z 



l + z /L 4 



- In z ln 2 (l - z) - 2 In z ln(l - z) ln(l + z) + - In 2 z ln(l - z) 
-2 ln(l - z)Li 2 (l - z) - 21n(l + z)Li 2 (l - z) + lnzLi 2 (l - z) 

- In zLi 2 (-z) + 2Li 3 (l - z) + S 12 (l - z) + S 12 (z 2 )])] . (A3) 

Here the factor (—1)™ originates from the non-planar diagrams (namely 13, 17 and 18 in 
figure 2 of 0). It multiplies that part of the matrix element which is needed for the mass 
factorization of physical processes with two identical quarks in the final state. The unphysical 
part (see Eq. ( |2.19| )) is equal to 
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NS,EOM 
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E' 



F C f 



4z - 2(1 - |) + e[2z - (1 - £)] [-1 + lnz + ln(l - z)] 
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_ „ r 20 88 . _ r 16 32 
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3 

16 1 16 
+— (- -3z-2z 2 )[lnzln(l + z) +Li 2 (-z)] +ln 2 z(12z + — z 2 ) 

O Z o 

+ 16Li 2 (l - z)(l - 3z) + 8z(l - z){-2 ln(l - z)((2) + 2 lnz((2) 
-In 2 zln(l - z) - 41nzLi 2 (l - z) - 21n(l - z)Li 2 (l - z) 
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+6C(3) + 6Li 3 (l - z) - 6S 12 (1 - z)}] 

8fi 4Qfi 4 1 fi 

+C A C F [-- + —z + 4C(2)(1 - 2z + -z 2 ) - - ln(l - z)(l + 4z) 

+ In z(6 - 28z) - \z 2 In 2 z - 4[ln z ln(l - 2) + 2Li 2 (l - 2)] (1 - 2z) 
8 1 

--(- -3z-2z 2 )[lnzln(l + z) + Li 2 (-z)] + 8z(l -z){-C(2)lnz 

- In 2 zln(l - z) + C(2) ln(l - z) + ln(l - z)Li 2 (l - 2) 
+21nzLi 2 (l - z) - 3C(3) - 3Li 3 (l - z) + 3S 12 (1 - z)}] 
+^7)^(5 - 16z) - jj(l - 2z){lnz + 2 ln(l - z)}] 



-y (-l) n (C| - - z + z 2 [-2C(2) + In 2 z\ 

- In z(l - 2z) + (~ - 3z - 2z 2 ) [In z ln(l + z) + Li 2 (-z)] } 

Now we give the pure singlet terms (see Eq. ( |2.21D ) 



.4 



PS.PHYS / 
99 



(z, , -) = F 2 n/CW>[^{5(l - z) + 2 In 2(1 + z)} 



rl6 



g 

+-{(1 - z)[7 + 10 ln(l - z)] + lnz(ll + 3z) 

+4(1 + z) [In z ln(l - z) + - In 2 z + Li 2 (l - z)] } 
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The result for Eq. ( ^22]) is 
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The OME in Eq. ( gjg ) is equal to 
-p 2 1 
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+ (l-2z){lnz + ln(l-z)} 2 



(A7) 



The next OME A gq can be split into physical and unphysical parts according to Eq. ( 2.25 ). 
The former becomes (see Eq. ( |2.26|) ) 
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The unphysical part given by Eq. (l2T27f) is equal to 
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-n 2 T 2 f^(l-2)} 

+^[-^(565 - 4032) + 2 -f-^- - C (2)(^ - fz + 82 2 

19 N x r 11141 214 w . 854 w . w 
-) + 5(1 - 2){— -C(2) - — C(3) + 5C(2) 2 } 



l-z' v /L 54 3 ^ v ' 9 

4 10 

+C(3)(-14 + 522 + 152 2 - — — - ) 

1 + 2 1 — 2 

+ ln(l - 2)C(2)(10 - 282 - 52 2 + -^-) + ^ ln 3 (l - 

L Z o 

x (1 - 22 + — !— ) + 2 ln 2 (l - 2)(-59 + 422 + -^-) 

1 — 2 1 — 2 
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16 2 

- ln(l - z)Li 2 (l - z)(7A + 108z + 5z 2 + ) 

1 + z 1 — z 

1 506 

-- ln(l - *)(-1267 + 959,2 + ) 

9 1 — z 

-4(5 + 8z + —— ){ln(l + z)£(2) + 2 ln(l + z)U 2 (-z)} 
1 + z 

-ln<(2)(30 - 12z - 5z 2 + -^—) - 41nzln 2 (l - z)(5 + 22z 

2 5 1 9T 

+- ) + -lnzln(l - z) (-245 - 29,2 + ) 

L -\- z L — z 6 L — z 

-4 In ^ In 2 (1 + z) (5 + 8,2 + - 2 In 2 * ln(l + z) (1 + 4z 

+ _l_ ) _ ln ^ ln(1 _, )(37+5 4 2 + ^- I l-- I 4-) 

_2 In 3 2 (12 - — + — — ) - In 2 z(3l + Alz - Az 2 - —) 
1 + z 3 1 — z 1 — z 

-21nzLi 2 (l -z)(10 + 60^ + 5z 2 - - ^-) 

2 2 226 

-41nzLi 2 (-z)(l + 4^ + + ) - lnz(— + 74,2 

i + z 1 — z 3 

521 1 . . w 13 1 , 

+ __)-2L l2 (l-,)(-21 + 45, + T — ) 

+4(8 -- + 7z- 22 2 ){ln^ln(l + z) + Li 2 (-z)} 

z 

16 6 

+Li 3 (l - s)(62 + 132* + 15,2 2 + ) 

1 + z 1 — z 

+4Li 3 (-z)(l + Az + ^— ) + S 12 (l - *)(46 - 148,2 - 15,2 2 

+ ih + Th ~ ss " ( - :){7 + 12z + TT7> + 8(1 + 2z + TT7> 

x{-21n(l + z)Li 2 (l -z) -21nzln(l-z) ln(l + z) + Si 2 (z 2 )}] 

2 64 1 4 

+n / C A 7>[--(473 - 563,2) - — + -C(2)(13 + Wz + Az 2 

y o i — ^ o 

6 , r/i xr 3224 236 l36 >/oX , 
__) +i (l- 2 ) { __ + _ C (2) + _C(3)} 

-8 ln 2 (l -z)(l-2z+ — !— ) + £ ln(l - 2) (97 - I31z 

1 — z 9 

34 4 8 
+- ) - - lnzln(l - z)(ll - 22,2 + ) 

X -2 o X ,2 



2 2 4 37 1 
-4 In 2 *(3 - 5z + -z 2 + ) + - lnz(10 - 37^ + ) 

3 1 — z 3 31 — 2 

-f Li 2 (l - z)(3 - -?—) + ^(3 + - + 3z + ^ 2 ){ln^ln(l + z) 
3 1 — z o z 

+U 2 (-z)} + 8,2 2 {ln(l - z)U 2 (l - z) - In z((2) + ]- In 2 * ln(l - 2) 
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+ ln(l -z)£(2) + 21nzLi 2 (l - z) - 3C(3) - 3Li 3 (l - z) +3S 12 (1 - z)}] 
+n / C F T / [5(l - z){-- + 16C(3)} + — (1 -z)- 4C(2)(3 + 5z + -z 2 ) 
+40 In 2 (1 - z) (1 - z) - 88 ln(l - z) (1 - 2) + 8 In z ln(l - z) (3 - 7s) 



+2 In 2 z(7 -llz + -z 2 ) - - In z(l - 65z) 
3 3 



32, 1 n 

— 3 + - + 3z + z 2 

3 z 



x{lnzln(l + z) + Li 2 (-z)} + 16(1 + z)(21n(l - z)Li 2 (l - 2) 

1 3 7 
+-ln<(2) + lnzln 2 (l - z) + -ln 2 zln(l - z) H In 3 z 

+21nzLi 2 (l - z) - Li 2 (l - z) - 2Li 3 (l - 2) + S 12 (l - z))] 
+n 2 T f 2 165(l-^)[l-^C(2)]" • 



9 



(A10) 



Finally we have determined the unrenormalized Z^ s from Eq. ( 2.15|) . The numerator, 
which is equal to the spin averaged OME A^ 9 S,PHYS , is given in Eq. (A. 3) of ]E|. The 
denominator has been computed above (see Eq. (|A3[)). The result is equal to 



75,NS 



5(1 — z) + (i s S £ Cf 
+2z-(l-£)}" 
+a 2 S 2 [C 2 F {- 16(1 

+CUCU^(1-; 



5(1 - z) +e{ -4(1 - z)(ln(l - z) +lnz) 



8(2 + z) lnz + 16(1 - z) lnzln(l - z)\ 



562 460 , 

— + — z + 8(l-z)C(2) 

80 



+ y (1 - z) ln(l - z) + (| - yz) lnz - 4(1 - z) In 2 z} 
+n/^l>{ - if (1 - z) + f - f z - f (1 - z) ln(l - z) 
-^(l-z)lnz} 

+ (-l) n (C| - ic A C F )f 8(1 + z) (4Li 2 (-z) + 4 In z ln(l + z) 



+2C(2) -ln 2 z-31nz) -56(1-2)} 



(All) 
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